On a relation of pseudoanalytic function theory 
to the two-dimensional stationary Schrodinger 
equation and Taylor series in formal powers for 

its solutions 

Vladislav V. Kravchenko 
Seccion de Posgrado e Investigacion 
Escuela Superior de Ingenieria Mecanica y Electrica 
Institute Politecnico Nacional 
C.P.07738 Mexico D.F., 
MEXICO 
phone: 525557296000, ext. 54782 
e-mail: vkravchenko@ipn.mx 

February 7, 2008 



Abstract 

We consider the real stationary two-dimensional Schrodinger equation. 
With the aid of any its particular solution we construct a Vekua equation 
possessing the following special property. The real parts of its solutions 
are solutions of the original Schrodinger equation and the imaginary parts 
are solutions of an associated Schrodinger equation with a potential having 
the form of a potential obtained after the Darboux transformation. Using 
L. Bers theory of Taylor series for pseudoanalytic functions we obtain a 
locally complete system of solutions of the original Schrodinger equation 
which can be constructed explicitly for an ample class of Schrodinger equa- 
tions. For example it is possible, when the potential is a function of one 
cartesian, spherical, parabolic or elliptic variable. We give some examples 
of application of the proposed procedure for obtaining a locally complete 
system of solutions of the Schrodinger equation. The procedure is algo- 
rithmically simple and can be implemented with the aid of a computer 
system of symbolic or numerical calculation. 

Keywords: Schrodinger equation, factorization, pseudoanalytic func- 
tions, generalized analytic functions, p-analytic functions, Darboux trans- 
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1 Introduction 



The appearance of a mathematical theory depends a lot on the preferences and 
mathematical tastes of its creator. Different people regard and describe the 
same mathematical results and ideas from different viewpoints, and this is one 
of the somces of the richness and progress of our science. 

The pseudoanalytic function theory is one of the clearest confirmations of 
this assertion. It was independently developed by two prominent mathemati- 
cians I. N. Vekua and L. Bers with coauthors, and presented in their books 
and |12| . The theory received further development in hundreds of posterior 
works (see, e.g., the reviews j2j and and historically it became one of the 

important impulses for developing the general theory of elliptic systems. Here 
the Vekua theory played a more important role due to its tendency to a more 
general, operational approach. L. Bers tried to follow more closely the ideas of 
classical complex analysis and paid more attention to the efficient construction 
of solutions. Among other results L. Bers obtained analogues of the Taylor 
series for psedoanalytic functions and some recursion formulas for constructing 
generalizations of the base system 1, z, z^,. . . . The formulas require knowledge 
of the Bers generating pair (two special solutions) of the corresponding Vekua 
equation describing pseudoanalytic functions as well as generating pairs for an 
infinite sequence of Vekua equations related to the original one. The necessity 
to count with an infinite number of exact solutions of different Vekua equations 
resulted to be an important obstacle for efficient construction of Taylor series 
for pseudoanalytic functions. 

Nevertheless as we show in the present work the Bers recursion formulas 
seem to have been specially designed for obtaining Taylor-type series for solu- 
tions of two-dimensional stationary Schrodinger equations admitting particular 
solutions which enjoy a peculiar property called in this work Condition S. The 
class of such Schrodinger equations is really wide. We show that if the poten- 
tial in the Schrodinger equation is spherically symmetric or it is a function of 
one cartesian, or parabolic, or elliptic variable, the corresponding Schrodinger 
equation belongs to this class. Moreover, the above mentioned cases are only 
some few examples. In general the Schrodinger equation belonging to the class 
must not necessarily admit separation of variables. 

The main result of the present work is a relatively simple procedure which 
allows us to construct explicitly a locally complete system of solutions of the 
Schrodinger equation by one known particular solution. Here the local complete- 
ness is understood in the sense that any solution of the Schrodinger equation 
can be approximated arbitrarily closely by a linear combination of functions 
from this system in a neighborhood of any point of the domain of interest. The 
global completeness is an open question nevertheless the possibility to obtain 
a sequence of exact solutions with such a special property as the local com- 
pleteness can be useful in different applications including qualitative analysis 
of solutions and numerical solution of boundary value problems. The main re- 
sult is based on a chain of observations some of them representing independent 
interest. 
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First of all we observe (Subsection 3.1) that given a particular solution of the 
stationary two-dimensional Schrodinger equation the corresponding Schrodinger 
operator can be factorized just as in a one-dimensional situation. In the consid- 
ered two-dimensional case the factorizing terms are operators + ^J^^ ^^'^ 
dz — ^j^C, where fo is a particular solution of the Schrodinger equation 

(-A + i.)/ = (1) 

and C is the complex conjugation operator. This observation gives us a simple 
relation between the Schrodinger equation and the Vekua equation 

(h+^c)w^O. (2) 
Jo / 

Every solution of one of these equations can be transformed into a solution of 
the other and vice versa. 

Next we show that solutions of this equation are closely related to solutions 
of another Vekua equation 

d^-^c)w ^0 (3) 

Jo / 

which we call the main Vekua equation. For this equation we always have a 
generating pair {F, G) in explicit form and the {F, G)-derivative of W (the op- 
eration introduced by L. Bers) is a solution of (01 ■ The {F, G)-antiderivative of 
w is a solution of Q. Moreover, the real part of W is necessarily a solution of ^ 
and the imaginary part of is a solution of another Schrodinger equation with 

the potential + ^ ^ which is precisely the potential which would 

be expected to obtain after the Darboux transformation (see, e.g., [?])■ We 
obtain a transformation which allows us to construct the imaginary part of W 
by its real part and vice versa obtaining in this way an analogue of the Darboux 
transformation for the two-dimensional Schrodinger equation. Here we should 
say that this transformation is not yet a long-sought definitive solution of the 
problem of a multidimensional generalization of the one-dimensional Darboux 
transformation (see, e.g., ^U]) because it is not clear how to include in our con- 
sideration the eigenvalues of the operator. Nevertheless it is a certain progress 
in generalizing the Darboux transformation and deserves more attention. 

Let /o be a function of some variable p ■ fo = fo{p) such that the expression 
Ap/ |Vp|^ is a function of p. We denote it by s{p) = f^p ^'^'^ thai fo 
satisfies Condition S. We show that under this condition for equation Q the 
formal powers in the sense of L. Bers can be constructed explicitly. It should be 
noted that formal powers play a part analogous to that of powers of the inde- 
pendent variable z in classical analytic function theory. They give us analogues 
of Taylor series expansions for solutions of (PJ, and locally any solution of ^ 
can be approximated arbitrarily closely by a finite number of first members of 
its Taylor series. 
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As was explained above, the real parts of formal powers are solutions of the 
Schrodinger equation (Q) and therefore similarly to real parts of powers of z 
which are very important in theory of harmonic functions and are widely used 
for numerical solution of boundary value problems for the Laplace equation, the 
real parts of formal powers give us a locally complete system of solutions. 

For the sake of simplicity we consider the Schrodinger equation with a real 
valued potential, and in the last section we explain how our results can be 
generalized to the case of a complex valued potential. 



2 Some definitions and results from pseudoana- 
lytic function theory 

This section is based on the results presented in |3] and @]. Let $7 be a domain 
in R^. Throughout the whole paper we suppose that is a simply connected 
domain. 

2.1 Generating pairs, derivative and antiderivative 

Definition 1 A 'pair of complex functions F and G possessing in U, partial 
derivatives with respect to the real variables x and y is said to be a generating 
pair if it satisfies the inequality 

lm(FG) > infl. 

Denote dz ^ ^ + *^ ^^"^ ^ ^ ~ (usually these operators are intro- 
duced with the factor 1/2, nevertheless here it is somewhat more convenient to 
consider them without it). The following expressions are known as characteristic 
coefficients of the pair {F, G) 

FG^ - F^G , FG^ - F^G 



FG~FG ' ' FG~FG 



FGz ~ F-xG ^ FGz ~ FtG 



A(F,G) = 7= =7r' B^F,G) 



FG-FG ' ' ' FG-FG 

where the subindex z or z means the application of or dz respectively. 

Every complex function W defined in a subdomain of Q admits the unique 
representation W = (j)F + ipG where the functions (j) and V' are real valued. 
Sometimes it is convenient to associate with the function W the function lu = 
4> + i^. The correspondence between W and lo is one-to-one. 

The {F, G')-derivative W — '^'^^'^'^ of a function W exists and has the form 
W = (f>zF + ^PzG^Wz~A^F,G)W-B^p,G)W (4) 
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if and only if 

chF + V^-G = 0. (5) 
This last equation can be rewritten in the following form 

= a{F,G)W + 6(F,G)W^ 

which we call the Vekua equation. Solutions of this equation are called (F, G)- 
pseudoanalytic functions. liW is (F, G)-pseudoanalytic, the associated function 
LO is called (F, G)-pseudoanalytic of second kind. 

Remark 2 The functions F and G are (F, G)-pseudoanalytic, and F = G = 0. 

Definition 3 Let (F, G) and (Fi,Gi) - be two generating pairs in fl. (Fi,Gi) 
is called successor of (F, G) and (F, G) is called predecessor of (Fi, Gi) if 

a(Fi,Gi) = 0-iF,G) O-nd &(Fi,Gi) = -B{F,G)- 

The importance of this definition becomes obvious from the following state- 
ment. 

Theorem 4 Let W be an {F,G)-pseudoanalytic function and let (Fi,Gi) be a 
successor of (F, G). Then W is an {Fi,Gi)-pseudoanalytic function. 

Definition 5 Let{F,G) be a generating pair. Its adjoint generating pair {F, G)* • 
(F*,G*) is defined by the formulas 

2F 2G 
F = = — ^— , G = 



FG-FG FG-FG 
The (F, G)-integral is defined as follows 

j Wd(F,G)Z^^[F{zi)Ke j G*Wdz + G{zi)Yie j F*Wdz 

where r is a rcctifiablc curve leading from zq to zi . 

liW = (j)F + ^G is an (F, G)-pseudoanalytic function where and ^ are 
real valued functions then 

/ Wd(F,G)Z = W{z)-(l){zo)F{z)-i,{zo)G{z), (6) 



and as F =G = 0, this integral is path-independent and represents the (F, G)- 
antiderivative of W. 
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2.2 Generating sequences and Taylor series in formal pow- 
ers 

Definition 6 A sequence of generating pairs {{Fm, Gm)}, m = 0, ±1, ±2, . . . , 
is called a generating sequence if {Fm+i,Gm+i) is a successor of {Fjn,Gm)- If 
(Fo,Go) = {F,G), we say that {F,G) is embedded in {{Fm,Gm)}- 

Theorem 7 Let {F, G) he a generating pair in f2. Let fli he a hounded domain, 
Oi c O. Then {F,G) can he embedded in a generating sequence in fli. 

Definition 8 A generating sequence {{Fm, Gm)} is said to have period /j. > if 
{Fm+ii,Gm+ij,) is equivalent to {Fm,Gm) that is their characteristic coefficients 

coincide. 

Let W be an (F, G)-pseudoanalytic function. Using a generating sequence 
in wliicli (F, G) is embedded we can define tlie liigher derivatives of W by the 
recursion formula 

W^'^=W; = ferf^!!, m = l,2,.... 

dz 

Definition 9 The formal power (a, zo; z) with center at zo S ft, coefficient 
a and exponent is defined as the linear combination of the generators Fm, Gm 
with real constant coefficients X, /x chosen so that XFm{zo) + iJ,Gm{zo) = a. The 
formal powers with exponents n = 1,2, . . . are defined by the recursion formula 

Z^^+'\a,zo;z) = (n+ 1) T ^^"[^(a, zq; C)ti(F„,G^)C- (7) 

This definition implies the following properties. 

1. Zm\a, Zo; z) is an {Fm, G„i)-pseudoanalytic function of z. 

2. If a' and a" are real constants, then Zm\a' +ia" , Zq; z) = a' Zm\l, Zq; z) + 
a"Z^\i, zo;z). 

3. The formal powers satisfy the differential relations 

^ '-^ =nZ^^^^'{a,zo;z). 

4. The asymptotic formulas 

Z^\a,zo;z) ^ a{z- Zo)'', z ^ zo 

hold. 
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Assume now that 

oo 

M^(z) = ^Z(")(a,zo;^) (8) 

where the absence of the subindex m means that ah the formal powers corre- 
spond to the same generating pair {F,G), and the series converges uniformly 
in some neighborhood of zq. It can be shown that the uniform limit of pseu- 
doanalytic functions is pseudoanalytic, and that a uniformly convergent series 
of {F, G)-pseudoanalytic functions can be {F, G)-differentiated term by term. 
Hence the function in |(SJ) is {F, G)-pseudoanalytic and its rth derivative 
admits the expansion 

oo 

W^"\z) = ^ n(n - 1) • • • (n - r + l)Z("-'-)(a„, zq; z)- 



From this the Taylor formulas for the coefficients are obtained 

W^["l(zo) 



(9) 



Definition 10 Let W{z) be a given (F, G) -pseudoanalytic function defined for 
small values of \z — zq\. The series 

oo 

5]z(")(a,zo;z) (10) 

with the coefficients given by |3) is called the Taylor series of W at zq, formed 
with formal powers. 

The Taylor series always represents the function asymptotically: 

N 

l¥(z)-^Z(")(a,Zo;^) = 0(1^-^01"^+'), z^zo, (11) 

ri=0 

for all N. This implies (since a pseudoanalytic function can not have a zero 
of arbitrarily high order without vanishing identically) that the sequence of 
derivatives {W^["l(zo)} determines the function W uniquely. 

If the series (|1U|) converges uniformly in a neighborhood of zo, it converges 
to the function W. 

Theorem 11 The formal Taylor expansion fl(J\) of a pseudoanalytic function 
in formal powers defined by a periodic generating sequence converges in some 
neighborhood of the center. 

We will use also the following observation from ^ p. 140]. 
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Proposition 12 Let b be a complex function such that bz is real valued, and 
let W = u + iv be a solution of the equation 

W- = bW. 

Then u is a solution of the equation 

dzdzu ~ (bb + bz)u ^ (12) 
and V is a solution of the equation 

&zdzV-{bb-bz)v^Q. (13) 



3 Relationship between generalized analytic func- 
tions and solutions of the Schrodinger equa- 
tion 

3.1 Factorization of the Schrodinger operator 

It is well known that if /o is a nonvanishing particular solution of the one- 
dimensional stationary Schrodinger equation 

(-^ + K.))/(.) = o 

then the Schrodinger operator can be factorized as follows 

dx^ \dx^foJ\dx foj- 

We start with a generalization of this result onto a two-dimensional situation. 
Consider the equation 

{-A + u)f = (14) 

in some domain C R , where A = ^2+^j2, and / are real valued functions. 
We assume that / is a twice continuously differentiable function. 
By C we denote the complex conjugation operator. 

Theorem 13 Let fo be a nonvanishing in f2 particular solution of il4)) - Then 
for any real valued continuously twice differentiable function ip the following 
equality holds 

{A-.)^^[9,+ ^C^(dz-^y. (15) 
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Proof. Consider 



^ Aip- ^p^ip = (A - I/) ip. 
Jo 



Remark 14 As ip in is a real valued function, we can add the conjugation 
operator in the second first-order operator on the right-hand side, and then M[^) 
takes the form 

The operator dz — ^jf^^; where / is the identity operator, can be represented 
in the form 

Jo 

Let us introduce the foUowing notation P = fadzfo^I- Due to Theorem 1131 
if fo is a nonvanishing solution of p4|l . the operator P transforms solutions of 
p4|) into solutions of the equation 

dz + ^c)w^O. (16) 
JO / 

Note that the operator dz applied to a real valued function ip can be regarded 
as a kind of gradient, and if we know that = $ in a whole complex plane 
or in a convex domain, where $ = <I>i + i$2 is a given complex valued function 
such that its real part $i and imaginary part <i>2 satisfy the equation 

9y$i+a,$2-0, (17) 

then we can reconstruct ip up to an arbitrary real constant c in the following 
way 

ip{x,y)^ $i(77,y)d??- / $2(2:0, O^^C + c (18) 
J^o Jyo 

where {xq, yo) is an arbitrary fixed point in the domain of interest. 
By A we denote the integral operator in (|18|) : 



A['P]{x,y) 



/ $i(?7,y)dry- / $2(a;o,Orf^ + c 
•Ixo Jyo 



Note that formula (|18(l can be easily extended to any simply connected domain 
by considering the integral along an arbitrary rectifiable curve F leading from 
{xo,yo) to {x,y) 



ip[x, y) = j ^idx - $2^2/ 
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Thus if $ satisfies (|17ll . there exists a family of real valued functions ip such 
that dz^ = given by the formula (p = A[^]. 

In a similar way we define the operator A corresponding to dz- 

l-x ry 

Jxo Jyo 

Consider the operator 5 = foAfQ^I. It is clear that PS = I. 

Proposition 15 '6! Let fo be a nonvanishing particular solution of jJ^I ) and w 
be a solution of \lb\l . Then the function f = Sw is a solution of il4V - 

Proposition 16 '6/ Let f be a solution of \14}) - Then 

SPf = f + c/o 

where c is an arbitrary real constant. 

Theorem ^| together with Proposition E| show us that equation H14|l is 
equivalent to the Vekua equation 116|l in the following sense. Every solution of 
one of these equations can be transformed into a solution of the other equation 
and vice versa. 

3.2 The main Vekua equation 

Equation l|16|l is closely related to the following Vekua equation 

dz-^c)w ^0. (19) 



fo 

To see this let us observe that the pair of functions 

F^fo and G=-^ (20) 
Jo 

is a generating pair for H19|) . Then the corresponding characteristic coefficients 
^{F.G) ^-nd B(^p Q) have the form 

A OR ^^-^0 

Jo 

and the {F, G)-derivative according to is defined as follows 

JO V fo J 

Comparing B(^p Q) with the coefficient in (I16|l and due to Theorem 0] we 
obtain the following statement. 
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Proposition 17 Let W be a solution of 1^19^) . Then its {F,G)- derivative, the 
function w = W is a solution of i|_?6]] . 

This result can be verified also by a direct substitution. 
According to 10 and taking into account that 

i^*--^/o and G* - l//o, 

the {F, G')-antiderivative has the form 

w{Od(p,G)C - \ Re f - jfr Re ifoiC)w{OdC 

2 V Jzo /o(C) .fo{z) Jzo 



(21) 



and we obtain the following statement. 

Proposition 18 Let w be a solution of \lb\) . Then the function 



W{z) = / u;(C)rf(i.,G)C 



IS a s 



olution of 1^19\) . 



Let (j) and "0 be real valued functions. It is easy to see that the function 
= 4>fo fo is a solution of (|19|) if and only if (j) and ip satisfy the equation 
ip- — ifQ(f>z = which is equivalent to the system 

"fpx + /o 0y = 0, tpy ~ fl(j)x = 

defining so called p- analytic functions (see [H] and [HI) with p = f^. 

Proposition 19 Let W be a solution of Then u — ReW is a solution of 

{14)) and V = Im W is a solution of the equation 

^+„_2(!^)'Uo. (22) 

Proof. Observe that the coefficient b — in p9|l satisfies the condition of 
Proposition El 



bz = 



fa \ fo J \ fo 

Thus, according to Proposition u is a solution of H12() and u is a solution of 

l(T^ . Calculating the expressions bb + bz = v and bb — bz = 2 ^^^^ 
finish the proof. ■ 



2 

— ly we 
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Proposition 20 Let u be a solution of j Then the function 

veker |^A + ^-2(^^ 

such that W — u + iv is a solution of 03), is constructed according to the 
formula 

V = f,'A{tf^ch{f„-'u)). (23) 
It is unique up to an additive term c/q"^ where c is an arbitrary real constant. 



Given v € ker 




the corresponding u g ker (A — v) can 



be constructed as follows 

u = -faA{tf^^ck{fov)) (24) 

up to an additive term c/q. 

Proof. Consider equation l(T{?|l . Let W — (f)fo + itp/fo be its solution. Then 
the equation 

i^-ifo^^O (25) 
is vahd. Note that ii u ^ KeW then (f> = u//o- Given (p, ip is easily found from 

= Aizf^ch). 

It can be verified that the expression A{ifQ(j>z) makes sense, that is dx{f^4>x) + 

dyifl^y) =0. 

By Proposition ^1 the function v = /q ^-0 is a solution of (|22|l . Thus we 
obtain H23|l . Let us notice that as the operator A reconstructs the scalar function 
up to an arbitrary real constant, the function v in the formula H23|l is uniquely 
determined up to an additive term cf^^ where c is an arbitrary real constant. 

Equation H24|l is proved in a similar way. ■ 

Remark 21 The potential in the Schrddinger equation l^2S^) has the form of 
a potential obtained after the Darboux transformation (cf. and thus 

formulas f2y\) and can be considered as a two-dimensional analogue of the 
Darboux transformation, though it is not clear how to include in our consider- 
ation the eigenvalues of the operator. 

Remark 22 When v = Q and fo = 1, equalities \2!^) and \24^ turn into the 
well known in complex analysis formulas for constructing conjugate harmonic 
functions. 

Remark 23 Equation can be written as follows 

{foMo'P^ + ^I^^^hP-) W - (26) 
where P+ = I {I + C) and p- = ~ C). 
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The form of the operator in (|26|l suggests the following form of an inverse 
operator 

H^^\ (/o3 + 1/0-13 (z/o<&)) (27) 



where $ must be such function that the expressions ^ (/g I'l') and yl(i/o$) 
make sense. 

Proposition 24 The function W — defined by is a solution of l^2b\l 
and equivalently of fiy\) if and only if ^ is a solution of fj6|) . 

Proof. Assume that $ is such that the expressions A (fo^^) and A (i/o<&) 
make sense, that is 

dy Re - d., Im = (28) 

and 

dy Re (^/o$) - a, Im (z/o$) = 0. (29) 
In this case let us substitute the function W — in (|26|l . We have 

I {fodzfo'P^ + ^fo'^foP-) {foA (/g-i*) + */„-i3(z/o$)) 

= i (/oa^A (/o-i$) + z/-ia^A(*/o$)) = 1 ($-$) = 0. 

Now let us prove that (|28|l and (|29|l are equivalent to the fact that $ is a 
solution of (^BJ. Denote ^ = Re$ and ip = lm^. Then ||211l and can be 

written as follows 

and 

(/o0) + dy ifo^) = 0. 
The last two equalities are equivalent to the system 

d^(j) + dyip = — 

Jo Jo 

Jo Jo 

which can be rewritten as the equation 

/o 
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Remark 25 This proposition shows us that in the case of equation l^iyfl the 
{F,G)-antiderivative can be calculated using \21^ . Indeed, let W be a solution 

of ^lyfl - Consider its {F,G)- derivative W which is a solution of l^lb]} due to 
Proposition |j7| It can be written as follows 

w^fod4f^\)+if,'dAfov), 

where u — Re W and v — Im W . Consider 
1 .-1. 



HCW = - l^foA l^fo'CWj + l^ifoCW 
- lhA{d^{fo'u)) - foA{ifo'd^{fov)) 



+ tf.'A {ifl(h ifo'u)) + if^'A (d^ifov)) 
= \{u + tv^ foA (ifo^d^ifov)) + tfo'A {tfld^ {fo'y))) 
Frome here, due to Provosition MiA we obtain 



HCW ^u + iv + cifo+'^, 

fo 

where Ci and Ci are arbitrary real constants. 

Thus, application of the operator HC to solutions of jJ6)] gives us exactly 
the same result as the {F,G)-antiderivative defined by \21}) . 



4 Taylor series in formal powers for pseudoana- 
lytic functions and solutions of the Schrodinger 
equation 

In this section we show how for a quite ample class of potentials by one known 
particular solution of (|14(l one can always construct an infinite sequence of its 
solutions possessing the property of local completeness. 

4.1 Condition S 

Lemma 26 Let ip be a nonvanishing analytic function. Then solutions of the 
equations 

Wy = bW (30) 

and 

vjy — —bw (31) 

are related in the following way. If W is a solution of ^30\) then w — ipW is a 
solution of ^31\). and if w is a solution of jj'ij) then W — w/(p is a solution of 

nrh) . 
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The proof of this statement is obvious. 

Proposition 27 Let /o be a function of some real variable p: fo = fa{p) such 
that for some real valued nonvanishing function r/ the equation 

d^{rid,p) = mn (32) 

holds. Denote Lp = irjpz- Then 

dzfo ^ _<£_dzfo 
fo ^ /o ' 

and if W is a solution of ^19\} . the function w = ipW is a solution of ^16]) and 
vice versa, if w is a solution of H6]} . the function W — w/ip is a solution of 

Proof. Consider the expression = Observe that ippz ~ —fp-z- Then 

dzfo ^ _f_foPz_ ^ f dzfo 
fo ^ fo /o ' 

From (|32() it is evident that (p is analytic. Then denoting 6 ~ we see that 
H19|l is equation (|30|) from Lemma ESI and equation Hlt)|) is equation 1311) . Thus 
by Lemma EBl we obtain the result. ■ 

Consider the following condition introduced in 0. 

Condition 28 ( Condition S) Let fo be a function of some variable p ■ fo = 
fo{p) such that is a function of p. We denote it by s[p) = f^^- 

The following proposition gives us a description of all possible solutions of 

Proposition 29 For a real valued nontrivial function p there exists a real val- 
ued nonvanishing function rj such that |5'ij)) holds if and only if p satisfies Con- 
ditionW^ 

Proof. Let p satisfy Condition|2Hl Denote rj = e^'^ where S is the antiderivative 
of s with respect to p. Consider 

^ ivdzP) = dz (e-^pz) = -se-^ | Vp|' + e'^^Ap = 0. 



Now assume that for p there exists a real valued nonvanishing function ry 
such that (|32|l holds. Then Ap + ^Pz = or in another form: 

rjl _ Ap 

If p is harmonic then Condition 1281 is obviously fulfilled, so let us consider the 
opposite case assuming that p is not harmonic. The last equation can be written 
as follows 

V Inyy = /iVp, 
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where /i — — ■ In order that the product /iVp be a gradient it is necessary 

that [V/i X Vp] = which imphes that /i = ■ 

Thus, due to Proposition 1291 we have that the function ip from Proposition 
l^has the form ip = ier^ where S{p) — J dp. 

4.2 Some examples of functions satisfying Condition S 

Examples of variables p which satisfy Condition S are numerous and important 
in applications. As was mentioned above any harmonic function p fulfills the 
condition and obviously 77 = 1 in l|32|l . That is, for example, p{x, y) — aix + a2y, 
where ai and 02 are arbitrary real constants, or p(x, y) — xy belong to that class. 
An important example is p{x,y) = r — ^ x"^ +2/^. In this case = |- 

The parabolic coordinate p[x, y) = r + x also fulfills Condition S: 
Consider the elliptic coordinates p and 0: 

^ — 2 cosh p cos 9, y ~ 2 ^^'^^ ^ ^^^^ ^' 
It is convenient to consider the magnitudes 



ri = \/(a;+-^) +y2 = ^(cosh^ 



Qj\ Qj 

r2 = J(^a;--j + y^ ^ - (coah p ~ cosd) . 

Let us verify Condition S for instance for the variable p. It is somewhat easier 
to consider p = a cosh/i. If Condition S is fulfilled for p then it is obviously true 
for p. We have p = ri + 7'2, and 

z + a/2 z-a/2 

Pz ^ 



n r2 



Then 



z + a/2 z-al2\(z + a/2 z-a/2 



rir2 



and 



We obtain 



Ap. 



rir2 
Ap p 



\Vpf P' 

Thus in all considered cases when p is one of the cartesian coordinates, when 
p — r, when p is one of the parabolic coordinates or when p is one of the 
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elliptic coordinates Condition S is fulfilled. Moreover, as the Laplacian admits 
separation of variables in all the mentioned coordinate systems, then if the 
potential is a function of such p, there exists a particular solution of H14|l 
/o = fo{p): and the results of this section are applicable to all Schrodinger 
equations with potentials depending on such p. 

We should emphasize first that these are only some examples which definitely 
do not exhaust all interesting in applications situations that can be covered by 
Condition S. And second, in order that such a solution /o = /o(p) exist fulfilling 
Condition S, it is obviously not necessary that i/ be a function of p. 

4.3 Explicitely constructed generating sequence for the 
main Vekua equation with /o satisfying Condition S 

In what follows we assume that /o is a nonvanishing solution of l|14() satisfying 
Condition S. 

Theorem 30 Let ip — ie^^ pz in fl. Then the generating pair {F,G) 
with F = fo and G = i/fo is embedded in the generating sequence {F„i,Gm), 
m = 0, ±1, ±2, . . . with F,„ = ip'^F and Gm = V^G. 

Proof. First of all let us show that {F,n,Gm) is a generating pair for m = 
±1, ±2, .... Indeed we have 

Im(F„G™) = Im(^"F(p'"G) = Im(|(^|^" FG) > 0. 

Taking into account that (p- = it is easy to obtain the following equalities 

I 1 2m i~. 

«(i=^,„,G™) = m a(F,G) = 0, 

and 

-B(F,„_i,G,„-i) = :^^B(^F,G)- 

We should verify the equality 
which turns into the equality 

= hF,G) -Bi^F,G)- 

As 6(F,G) = ^jf^ and B(f^g) = by Proposition |53 we obtain that lO is 
true. Thus the sequence (-Fm, Gm), rn = 0, ±1, ±2, . . . satisfies the conditions of 
Definition 1^1 and therefore it is a generating sequence. ■ 

This theorem opens the way for explicit construction of formal powers of 
any order n > corresponding to the generating pair (/o,V/o) as well as to 
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any generating pair embedded in the sequence proposed in Theorem 1301 As a 
consequence L. Bers' theory of series expansion for pseudoanalytic functions can 
be used in order to obtain explicitly Taylor series in formal powers for solutions 
of the Schrodinger equation (fn|l due to Propositions [THI and HUi 

Due to we have that any pseudoanalytic function W can be approx- 
imated at any point zq € fl with an arbitrary precision by first N members 
of its Taylor series in formal powers. As any solution of H14|l is a real part of 
some pseudoanalytic function W satisfying it can also be approximated 

with arbitrary precision by the sum of real parts of the first members of the 
Taylor series in formal powers of the function W. 

Definition 31 Let u{z) be a given solution of \14{ l defined for small values of 
\z — zo\, and let W{z) be a solution of ^19\) constructed according to Proposition 
\^such that ReVF = u. The series 

oo 

^ReZ("'(an,zo;z) (34) 

with the coefficients given by 0j is called the Taylor series of u at zq, formed 
with formal powers. 

Theorem 32 

N 

u(z)-;^ReZ(")(a„,zo;;j) = o(|z-2or+') , z ^ zo, (35) 

ji=0 

for all N, and if the series J^j) converges uniformly in a neighborhood of zq, it 
converges to the function u. 

Proof, is a direct consequence of ■ 

Due to Theorem ll-jQI we are able to construct Z^"^^ (a„, zq; z) explicitly in many 
practically interesting cases. The simplest case is when /o is a function of one 
cartesian variable: /o = fo{y), that is we consider the equation 

-M{x,y) + y{y)f(x,y)^Q in a (36) 

Let us make a useful observation. In this case pz = —i and Lp = 1. Thus 
(Fm,Gm) = {F,G), m = 0,±1,±2,... , and obviously we have a periodic 
generating sequence with the period 1. Consequently, according to Theorem llll 
in the case under consideration we can guarantee not only the approximation 
Hll|l and H35|l but also the convergence of the Taylor series in formal powers in 
some neighborhood of the center. Let us consider the following example. 

Example 33 Let Q be the unit circle with center at the origin. Consider equa- 
tion iV6|) with v{y) = 6/ (y + 1) . One particular solution depending on y only 
can be chosen as follows fo{y) — (y + 1) . Lt is easy to find another solution 
of ^36]) . We choose it in the following form u{y) — {y + 1) . Next we con- 
struct the function v such that W = u -\- iv be a solution of Using \2!^) 
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we obtain v — 5x {y + 1) + c{y + 1) ^ where c is an arbitrary real constant. 
We choose it equal to zero, so v = bx (y + I) . It can be easily verified that 
the function W = {y + 1) + bix (y + 1) is indeed a solution of flU^) where 
— 3i{y+l)^^. Our aim is to find the Taylor series in formal powers of the 
function W at the origin. 

We find that W{0) = 1. Then using the definition of formal powers we find 

Z^^\l,0; z) — {y + 1)'^. In order to construct Z^^^ we need to calculate W{Q). 

We have W{z) = dzW{z) + Zi{y + l)-^W{z) = IQi {y + 1)"^ . Thus W{0) ^ m. 
In order to apply we find Z^^''{10i, 0; z) — IQi/ (y + I)''. Then using ^ and 
Remark \25\ we have 

(lOi, 0; z) = i/CZ(") {lOi, 0; z) = {y + ly^ - {y + if + 5ix {y + 

The function W{z) = Z^^^{1, 0; z) + Z^^\lOi, 0; z) should satisfy the equality 
W{z) — W{z) = 0(|z|^) when z ^ 0. We see that in fact W coincides with W as 

could be expected due to the fact that W is the generating function G multiplied 
by a real constant and hence W^"^^ = 0. 

Formal powers' Property 2 from Subsection 12.21 toerether with Theorem 1301 
allows us to construct in explicit form a locally complete system of solutions 
of (|19|l and of H14|) in the following sense. Due to Theorem 1201 when /q ful- 
fills Condition S we are able to construct the formal powers Z(")(l,zo;z) and 
zq] z), 71 = 0, 1, 2, . . . corresponding to any point zq G ^- Due to Property 
2 of formal powers we have that Z^'^\a, zq; z) for any Taylor coefficient a can be 
easily expressed through zq; z) and zq; z). Thus for any solution 

W of l|19|) there exists a linear combination of ^("'(1, zq; z) and Z'^"'\i, Z{)] z), 
n = 0, 1, 2, . . . iV such that (^TJ is valid. 

Hence for any solution u of (|14(l there exists a linear combination of Re zo; z) 

and ReZ(")(i, zq; z), ti = 0, 1, 2, . . . TV such that (jSU holds. 

Let us illustrate the procedure of construction of this locally complete system 
on the following example. 

Example 34 Let f2 be the unit circle with center at the origin and a, (3 positive 
real constants greater than 1. For the Schrodinger equation with 



we have the particular solution fo{x,y) = y^{x + a){y + P). Denote p ^ {x + 
a) {y + /3). Then Condition S is fulfilled and we have that the function ip from 
Theorem\^ is defined as follows (p = ie~^ pz = z + c, where c = a + if3. Let us 
construct the first formal powers Z(")(l,0;z) and Z^"\i,0;z). By Definition\^ 
we have 



Z(°)(l,0;z) = J (" + ")(f + ^^ and Z^^\^,Q, z) - ^ ^ "'^ 



a/3 V (a; + «)(?; + /?)■ 
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In order to construct Z^^\l,0; z) and Z^^"^ {i,0; z) by formula ^ we need first 
Z[°\l,0-,z) and z[°\i,0;z). We have that Fi ^ {z + c)y/{x + a){y + f3) and 
Gi = i{z + c)/y/{x + a){y + P). Then z["\l,0;z) = AiFi(z) + ^iGi(z) with 

Ai and /ii satisfying the equality Aii^i(O) + /iiGi(O) = 1 and z[^\i,0; z) = 
X2Fi{z) + fj,2Gi{z) with A2 and ^2 satisfying the equality A2i^i(0) + /i2Gi(0) ~ i. 
We obtain that 



and 



Consider 

Z«(a,0;z)= r zf\a,Q-Od(F.G)C 



\iV¥^a)WW) Re (a(c + c) + 



Im / ( A(C + c) (x' + a) (y' + /3) + m(C + c)*) dC,) 



y/[x + a){y + (3) Jo 

where A and fj, are real numbers such that A-F'i(O) + /iGi(O) = a and ( = x' + iy' . 
We have 

Re / (C + c)dC= / {{x^ - y'^)t + ax ~ I3y)dt ^ + ax - I3y, 

Jo Jo 2 

r (C + c)tdC ^ _ 2xyU^ayJj3^^^ 
^Jo {x' + a){y' + 13) Jo {xt + a){yt + (3) 

\a;(a?/ — /3a;) \ a / y{ay — (3x) \ (3 J J ay — [3x [3 [x + a) 

a(3 



= In 
and 



Im / (C + c)idC = Rc f (C + c)dC = ^ ^ + ax - Py, 
Jo Jo 2 

Im r {(:+c){x'+a){y'+f3)dC = ^-^ + ial3+xy)iay+px)+aPxy~- + ^^)' 
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The 



7(1). n ^ Vix + a){y + P) ,^, x'^~y'^ ^ f aP \ 

)(«, 0, z) = (A(^— + ax - (3y) + ^In (^^-—^^—^ j 



2 

+ m( h ax -/??/)). 



(A( — h (a/3 + xy){ay + /3x) + a/Jxj/ 



One can check that for any A and ^ the real part of this function is indeed a 
solution of jj^l ) with the potential J5'7| ). Substituting A and /i /rom or 
we obtain Z^^'^ (1, 0; z) anrf Z(^)(i,0;2:) respectively. 



5 Complex potentials 

Our approach can also be applied to the Schrodinger equation (|14|) with v being 
a complex function, though in this case complex numbers become insufficient, 
and one should consider the bicomplex generalization of the pseudoanalytic 
function theory. 

Together with the imaginary unit i let us consider another imaginary unit 
j such that i"^ = —1 and ij = ji. Bicomplex numbers have the form z\ + z^i 
where z\ and 22 can be considered as complex with respect to the unit j: zi,2 = 
2^1,2 + i2/i,2, 2:1,2 and y\^2 being real numbers. Bicomplex numbers obviously 
form a commutative algebra which contains a subset of zero divisors. 

Now we assume that v = v\-\- jv-i where and V2 are real valued functions. 
The factorization p5|l remains valid for (/? = + j'^i as well as all the results 
of the present work. 



6 Conclusions 

In the present work we considered the real stationary two-dimensional Schrodinger 
equation. With the aid of any its particular solution we construct the Vekua 
equation possessing the following special property. The real parts of its so- 
lutions are solutions of the original Schrodinger equation and the imaginary 
parts are solutions of an associated Schrodinger equation with a potential hav- 
ing the form of a potential obtained after the Darboux transformation. After 
having applied L. Bers' approach to this Vekua equation we obtained a locally 
complete system of solutions of the original Schrodinger equation which can be 
constructed explicitly for an ample class of Schrodinger equations, namely when 
the Schrodinger equation admits a particular solution satisfying the proposed 
Condition S. We established that in such special cases as of the potential being 
a function of one cartesian, spherical, parabolic or elliptic variable the condition 
is fulfilled. We gave some examples of application of the proposed procedure for 
obtaining a locally complete system of solutions of the Schrodinger equation. 
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The procedure is algorithmically simple and can be implemented with the aid 
of a computer system of symbolic or numerical calculation. The obtained sys- 
tem of solutions is a good candidate for numerical analysis of boundary value 
problems for the Schrodinger equation. 
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